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THEORETICAL  AND  PRACTICAL  ASPECTS  OF  SINGULARITY 

AND  EIGENMODE  EXPANSION  METHODS 

* 

A .  G .  Ramm 


5 1 .  Ir.f.  roducti  on . 

Vast  literature  has  been  written  on  singularity  a d  oigen- 
modc  expansion  methods  during  the  last  decade-  Engineers  and 
physicists  stimulated  interest  in  the  subject  (see  [1],  [2] , 

[3]  and  references  given  in  this  review).  Mathematical  analysis 
of  the  problems  was  initiated  in  [4],  [5]  and  was  pushed  con¬ 

siderably  further  by  M.  Agranovich  (see  [1]).  Nevertheless 
many  questions  in  the  theory  are  open  and  of  considerable  .in  to  in¬ 
to  engineers  and  mathematicians.  The  purpose  of  this  paper  can 
be  summarized  as  follows:  We  are  going  to  explain  in  a  simple 
way  the  principal  features  of  the  singularity  and  eiyennodo 
expansion  methods  and  to  formulate  explicitly  1)  what  has  been 
used  by  engineers  with  proof,  2)  what  is  important  for  practice 
3)  what  has  been  rigorously  established  and  ■',)  what  are  tlu- 
unsolved  mathematical  problems  in  this  field. 
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The  main  results  obtained  in  scalar  wave  scattering  theory 
were  generalized  to  electromagnetic  wave  scattering  without 
much  difficulty.  That  is  why  and  also  for  simplicity  we  restrict 
ourselves  by  the  presentation  of  the  theory  for  scalar  wave 
scattering. 


§2.  What  are  singularity  and  oiqenmode  expansion  method: 
(  SEM  and  EEM)  ? 

1.  What  is  EEM? 

Consider  the  problem 


(V2  +  k2)u  =  0  in  P  =  IR3  D,  k2  >  0  , 


(1) 


u|r  =  0 


(2) 


u  = 


u0  +  V 


—  exp{ik  (l.i  ,x)}  +  v 


(3) 


and  v  satisfies  the  radiation  condition 


<W 

a  t  >n~ 


-  i  k  v  =  o ( I x  )  as  I x I 


(4) 


Here  r  is  the  smooth  surface  of  a  finite  obstacle  D. 
If  we  look  for  a  solution  of  (1)  -  (4)  in  the  form. 


cxp(ikr  .) 
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f  ( t )  d  t ,  r 


x  t 


t|  . 


(5) 
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then 


I  ii.uwi  !  1 1  mu  ,  iu.  i  ii  ijw*' 1  IW«- 


A(k)f  5 


f  exp(ikrsfc) 


4  Trr 


st 


f(t)  dt  =  -  uQ(s)  ,  s  6  T.  (6) 


The  EEM  method  can  now  be  explained  as  follows.  Suppose  that 
the  set  of  eigenvectors  of  operator  A(k) 


A( k)$.  =  (k)  ,  j  =  1,2, .  (7) 

2 

forms  a  Riesz  basis  of  L  ( F)  =  H.  This  means  that  any  g  €  H 
can  be  expanded  in  the  series 


9 =  I  ' 

j=l  3  3 


(8) 


and 


ci  llgll2  i  l  |g-il2 

j=i  3 


-  C2 


c^  >  0 


(9) 


where  (|g||  is  the  norm  in  l  (f).  The  inequality  (9)  substi¬ 
tutes  the  Parseval  equality  for  orthonormal  bases.  A  complete 
system  in  H  does  not  necessarily  form  a  basis  of  II  (example: 

H  =  L2[0,1],  the  system  ^(x)  =  x-5 ,  j  =  0,1,2,...:  Not 
2 

every  g  e  L  (0,1]  can  be  expanded  in  the  series  c(x)  =  }j-Q  9j 

If  the  assumption  made  is  true,  then  equation  (6)  can  be 
solved  by  the  Picard  formula 

io  u 


f(s)  =  -  ) 


j=l  Xj(k) 


Sl__  .  (s) 
3 


(10) 
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where  the  coefficients  are  uniquely  defined  by  the  equality 
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u> 
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This  method  of  solution  of  the  scattering  problem  (i; - q 
is  called  EEM.  It  was  used  without  mathematical  analysis 
by  engineers  |1],  (3).  The  questions,  which  immediately 

arise,  can  be  formula  ted  as  follows:  1)  Is  it  true  that  the 
nonselfadjoint  operator  A(k)  has  eigenvectors  (e.g.  Volterra 
operator  has  no  eigenvectors) ;  2)  Is  it  true  that  the  set  of 
eigenvectors  of  A(k)  forms  a  Riesz  basis  of  H;  3)  suppose 
that  the  set  of  eigenvectors  (e  eigensystem)  of  A(k)  does 
not  form  a  Riesz  basis  of  H.  Is  it  true  that  the  root  system 
of  A(k)  forms  a  Ries2  basis  of  H  ? 

Let  us  explain  the  root  system.  Let  A  be  a  linear  operator 
on  H  A<p  =  A4>,  $  7*  0.  Consider  equation  A^  =  <j> .  If 

this  p":  c:;;.  is  solvable,  <)>^  is  called  a  root  vector  of  A 
corresponding  to  eigenvalue  A  and  eigenvector  4> .  If  exists 

consider  equations  A^  —  A k  >  1.  It  is  known  [7] 
that  only  a  f  j  n  i  to  number  r  of  root  vectors  4.^,. associated 
with  ([)  exist.  The  cha  i  »  ( 4,,  4^  , .  .  .  ,  4.^)  is  caLlod  a  Jordan  chain 
with  the  length  r  +  1.  The  union  of  all  root  vectors  of  a 
linear  operator  A  corresponding  to  all  eigenvalues  of  A  is 
called  the  root  system  of  A.  It  is  well  known  from  1 inear  algebra 
that  the  eigensystem  of  a  nonselfadjoint  operator  may  not 
form  a  basis.  For  example  if  the  operator  A  is  an  operator 


-4- 


then  A  has  only  one 


2  (  x 
xn  IR  with  the  matrix  A  =  I 

1 

eigenvector  so  that  the  eigensystem  of  A  does  not  form  a 
2 

basis  of  IR  .  It  is  also  known  that  the  root  system  of  any 
linear  operator  (matrix)  in  IRn  forms  a  basis  of  IRn  .  Of 
course  in  IRn  any  basis  is  Riesz  basis.  In  a  Hilbert  space 
(infinite  dimensional  space)  this  is  not  true.  For  practice 
it  is  important  to  have  affirmative  answer?  for  questions  1) 
and  3).  Indeed,  if  the  eigensystem  of  A(k)  does  not  form  a 
Riesz  basis  but  its  root  system  forms  a  Riesz  basis  of  H,  then 
it  is  still  possible  to  solve  equation  (6)  using  the  root 
system  of  A. 

2.  What  is  SEM? 

In  order  to  explain  what  SEM  is,  consider  the  problem 


u 


tt 


-  Au 

oo 


where  f  6  CQ  (a 
If  G(x,y,-p 

(-A  +  p2)G  « 

then  the  solution 


=  0  in  fl  ;  ujr=  0;  u|t=0  =  0,  ut|t=0 

smooth  function  which  vanishes  for  large 
2 

)  is  the  Green  function  of  the  problem 
5(x  -  y)  in  fi  ,  G]j-  =  0,  Rep  >  0 

of  (12)  can  be  written  as 

,c  +  i®  __ 

exp(pt) u(x,p)  dp, 

•  c  —  i«j 


=  f(x), 


(13) 


( 1  ' 


(12) 
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where 


u(x,p)  =  |  G(x,y,-p2) f (y)dy  (15) 

fi 

Suppose  that  u(x,p)  is  meromorphic  in  p  on  the  whole 
complex  plane  p  (this  is  actually  true  [8]- [12])  and  suppose 
the  following  estimate  is  valid 


u(x,p)  j  < 


1+  P 


,  a  >  1/2,  c  =  const,  >  0,  Rep  >  -A, 


(16) 


Imp  >  N 


A  ' 


where  A  >  0  is  arbitrary. 

This  estimate  fellows  from  the  Lax-Phillips  result  [12] 
and  the  arguments  given  in  [9],  [10],  [13],  Under  the  assump¬ 
tions  mace,  the  contour  of  integration  in  (14)  can  be  moved  to 
the  left  and  one  gets 


N  b.-l 

u  =  l  c^t  3  exp(p^t)  +  o ( exp ( -Re  p  t) ) 


j  =  l 


1 


1 


(17) 


where  p^  arc  the  poles  of  G(x,  y,  -p  ) ,  b^  are  the  mul¬ 
tiplicity  of  the  pole  p ^ .  Expansion  (17)  is  called  the  SEM 
expansion.  Actually  such  expansions  were  known  for  a  long  time 
for  various  concrete  problems  of  mathematical  physics  (especially 
in  cases  when  the  solution  can  be  represented  explicitly  in 
the  form  of  series).  The  main  difficulty  .is  to  prove  estimate 
(16)  which  allows  us  to  move  the  contour  on  integration  ; r  (14). 
If  only  the  mermorphic  nature  of  n(>:,p)  a  -  a  ‘ 


(17)  can 


is  established,  then  no  .....  a  ns  ion-  of  the  typo 

be  proved  in  yen era 1  because  there  is  a  possibility  that  poles 
-r  +  i„;  with  very  small  r  >  0  and  very  large  can 

exist..  The  general  Kit tag-LOf f  1  er  representation  can  not  be 
applied  for  derivation  of  SEN  expansion  (17),  because  this 
representation  uses  special  system,  of  expanding  contours,  while 
the  derivation  of  (17)  requires  the  possibility  to  move  our 
particular  contour  (c  -i<v,  c  +i»)  to  the  left.  From,  a  prac¬ 
tical  point  of  view,  the  SEM  expansion  is  used  at  present 
according  to  the  following  scheme: 

Suppose  that  only  a  few  terms  in  (17)  are  o&senf : ial,  e.g. 

1-3.  This  will  be  true  if  |  Re  p_.  |  »  Rep^  for  -j  >  3.  Then 
in  experiments  the  transient  field  u(x,t)  is  measured  a  ml 
each  p j ,  j  =  1,2,3  is  determined.  It  is  assumed  that  the 

p 

location  of  these  complex  poles  of  the  Green  function  G(x,y,-p‘) 
can  give  information  enough  to  identify  the  obstacle  (the 
scatterer  D) .  This  assumption  has  not  been  backed  theoretically. 
Nevertheless,  if  there  is  a  finite  set  of  scatters  (say  flying 
targets)  it  is  possible  to  believe  that  a  one  to  one  corre¬ 
spondence  can  be  established  empirically  between  the  scattorers 
and  the  correspond  i  ng  complex  pole:,. 

An  interest  1  r.g  in s:  pro:  !  -  ..i  can  be  formulated  in  con¬ 
nection  with  lb.- 

Inverse  problem::  Given  a  cr- 1  of  cor.  r  '  :  1  •<>  .  >  , 

I 

Ro  !j.  <  0 ,  i  -  'it  no..  i  '■  the 

'  J 

Green  fun.rt.ic:  cm.,-  .  ..  r-  •  s  w  .  i  h.  ••  •  a.'1.:- 

‘  P  •  ■■  7  bows  t . !  set  '  ■!.  i  .  .  i  v  -  ;■  1  ■  . ■  ; 


*.  ......  * 


Khat  restrictions  must  be  imposed  on  the  set  { p_^  )  ,  Re  <  0 
in  order  that  this  set  will  be  the  set  of  complex  poles  of  the 
Green  function  of  a  scatter'd*? 

If  the  scatterer  is  a  star-like  body  and  the  boundary  con¬ 
dition  is  the  Dirichlct  condition  then  the  set  { p . }  must 
satisfy  the  condition  |  Re  p ^  [  >  a  In  [  Iir.  p  |  +  b,  a  >  0  [12]. 
It  seems  that  no  other  information  on  the  problem  is  available. 

From  a  practical  point  of  view  this  problem  may  not  be  so 
important  as  it  seems.  First,  only  a  few  complex  poles  are 
available.  It  seems  hopeless  to  make  any  general  conclusions 
about  the  scaterer  from  this  information  without  severe  restric¬ 
tions  on  the  set  of  scatterers.  (For  example,  if  it  is  apriori 
known  that  the  scatterer  is  a  ball,  it  is  possible  to  determine 
its  radius  from  the  above  information.)  That  is  why  this  author 
thinks  that  from  a  practical  point  of  view  in  order  to  use  the 
SEM  for  identification  of  scatterers  it  is  more  useful  to  work 
out  tables  of  responses  of  the  typical  scatterers,  then  to  try 
to  develop  a  theory  of  the  posed  (which  is  very  interesting 
from  a  theoretical  point  of  view)  inverse  problem. 


with  the-  HEM  .!:.■■■ 


.ons  arise  naturally  in  connection 


•  integral  operators  in 
form  a  Riesz  basis  of  K? 


net  ci-..  with  th-. 


1) 


S 


roo  l 


3)  Do  the  complex  poles  of  the  Green  function  depend 


continuously  on  the  obstacle?  In  more  detail:  suppose  that 
Xj  =  Xj(t^,t2 ),  1  £  tl,t2  —  *  are  Parametri°  equations  of 
r,  yj  =  x^(t1#t2)  +  cZj(t1,t2),  0  £  t1,t2  <_  1,  0  <  e  <  1 

are  parametric  equations  of  the  surface  of  perturbed  scatterer. 

2 

Let  us  assume  that  x^  (t)  ,  z_.  (t)  6  C  (A),  t  =  t^t^  a  = 

{t^  £  t  ^ , t 2  £  1}.  Let  us  fix  an  arbitrary  number  R  >  0  and 
let  p j ,  1  £  j  £  r(R)  be  the  complex  poles  of  the  unperturbed 
Green  function  which  lie  in  the  circle  |p_.j  £  R-  Let  Pj(e) 
be  the  complex  poles  of  the  perturbed  Green  function.  Our 
question  can  now  be  formulated  as  follows:  is  it  true  that 
P(e)  ->  P j ,  as  e  -v  0,  uniformly  in  1  £  j  £  n(R)  provided 
that  the  numeration  of  p^(f.)  is  properly  done? 

4)  How  can  one  calculate  the  complex  poles? 

5)  What  are  sufficient  conditions  for  the  validity  of  SEM 
expansion  (17)  ? 

6)  Is  it  possible  to  calculate  complex  poles  via  calcula¬ 
tion  of  zeros  of  some  functions? 


§3.  What  has  been  rigorously  established  in  EEH  and  SEM  methods 

In  this  section  we  give  answers  to  questions  l)-6)  of 
Section  2.3.  No  proofs  will  be  given  but  the  results  obtained 
will  be  formulated  and  references  will  be  given.  Proofs  are 
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omitted  £ o c  thro**  rc-usor.s:  1)  they  arc  di f i  icul  t  fo:'  eft  usee rs, 
2)  they  are  long,  and  3)  they  can  found  in  the  cited  papers. 

1.  In  order  to  formulate  the  answer  to  question  1)  of  Section 

2.3  wo  must  explain  what  a  Ricsz  basi s  with  brackets  is.  Lot 

if.}  bo  an  or  then  or  mol  basis  of  H,  {h.f  be  a  cot.nl  etc  and 
'  J  3 

minimal  system  in  H.  (A  complete  system  {h^}  is  called  minir. 

if  the  system  {h.}"\h  is  not  complete  for  any  k,  k  -  1,2,3. 

3  k 

In  other  words  if  we  remove  any  element  of  our  system,  we 

obtain  an  incomplete  system)  .  Lot  <  rn^  <  .  .  -  be  an  arbi  ¬ 

trary  increasing  sequence  of  integers;  F..  is  the  linear  space 

with  the  basis  {f  ,  f  ,...,f  ,  H.  is  the  linear 

Ir,j-1  "j-1+1  “j"1 

space  with  the  basis  {h  , }  .  Suppose  that  there 

3~1  3 

exists  a  linear  bounded  operator  B  with  bounded  B  '  ,  defined 
on  all  H,  such  that  Bib  -  fb  ,  ;j  ~  1,2,...  .  Then  the  system 
{ h ^ }  is  called  a  Rios/.  basis  of  H  with  brackets.  This  defi¬ 
nition  is  equivalent  (sec  [14  1  )  to  the  following.  Let  P..  be 
projectors  in  II  onto  lb  .  Suppose  that  for  any  f  e  II, 

c,  llf2ll  i  l  li l*ic li2  1  c2  lif!!2'  c  >  0  , 

i  j=l  3 

then  the  system  [Jm]  is  called  a  Rios?,  basis  of  If  wii.h 
brackets . 

It  is  proved  that  the  root;  system  of  operator  A  ( k )  (see 
formula  (G)  )  forms  a  L  ies/  basis  v:i  •  h  brackets  ffij.  The  same 
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is  true  for  the  operator  arising  in  the  exterior  Neumann 
boundary  value  problem  (1] ,  [6] .  The  same  is  also  true  for  the 

electrodynamics  scattering  problem  [1J. 

2.  If  the  integral  operator  of  the  diffraction  problem  is 
normal  then  its  eigensystem  coincides  with  its  root  system. 

This  is  the  case  in  the  problem  ( 1 )  —  ( 4 )  for  a  spherical  surface 
T  and  for  a  linear  antenna.  First  this  was  observed  in  [5]. 

An  operator  is  called  normal  if  A*A  =  AA*,  where  A*  is  the 
adjoint  operator.  The  condition  AA*  =  A*A  is  the  condition  on 

T  provided  that  the  kernel  of  the  integral  operator  A  is  given, 
(see  [5]  for  details). 

3.  The  answer  to  question  3  from  section  2.3  is  affirmative 
(see  [6]  for  details). 

4.  A  general  method  (with  the  proof  of  its  convergence)  for 

calculating  the  complex  poles  of  the  Green’s  functions  in  dif¬ 
fraction  and  potential  scattering  theory  was  given  in  [4],  [15] 

(also  see  [6]).  The  method  can  be  explained  for  the  problem 
with  the  impedance  boundary  condition 

(V2  +  k2)u  =  0  in  ft  ,  -  h  u  |r  =  0  (18) 

where  h  =  const./  Re  h  >  0,  N  is  the  outer  normal  to  r  • 

We  look  for  a  solution  of  the  problem  without  sources  in  the 
form 
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- 1 
Jr 


dt 


From  (19)  and  (18)  it  follows  that 


o  =  Qa  , 


(19) 


(20) 


whore 


Qo  = 


'  _9_  exp(lkrst 

r  3Nt  27rrst 


-  a (s) ds  -  h 


exp (ikr  ) 


2  ir  r 


(s)cis  (21) 


St 


Let  {<5j}  be  a  Riesz  basis  of  H  =  L  (r), 


N 


°N  =  ' 


(22) 


Substituting  (22)  in  (20)  and  multiplying  by  in  H  we 

obtain : 

N 

Y  b .  .  (k)  c  .  =  0 ,  l<i<N;b..(k)=6---(Qi>.,  £■)•  (23) 

J  —  —  xj  ij  j  i 


This  system  has  nontrivial  solutions  if  and  only  if 


det  b^  (k)  =  0.  (24) 

The  left  -hand  side  of  (24)  is  an  entire  function  of  k.  Let 

be  its  roots.  Then  there  exist  the  limits  lim  k^  =  k 
rrt  „  ra  m 

N-x" 

and  k  are  the  poles  of  the  Green  function  corresponding  to 
m 

the  problem  ()8).  Moreover,  all  the  complex  poles  can  be  obtained 
by  this  method.  Proofs  are  given  in  [4],  (131,  [61.  From  a 
practice!  .  .  ;.>i:  of  view  there  are  two  nontrivial  points  in  por- 
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forming  this  method.  1)  caiculat i on  of  b ^ ( K )  by  formula. 
(23)  and  2)  numerical  solution  of  equation  (21).  For  bofh 
stops  there  are  methods  available  in  the  literature  on 
numerical  analysis. 


5.  Sufficient  conditions  for  the  validity  of  SEf!  expansion 
(17)  v.’ere  given  in  Section  2. 


6.  The  set  of  the  complex  poles  of  the  Green  function  of 
the  problem  (1)  -  (4)  coincides  with  the  set  of  the  complex- 
zeros  of  the  eigenvalues  p^(k)  of  the  operator  A(k): 


A(k)dn  =  Pn<kHn'  n  =  1,2,...  .  (25) 

Indeed,  let  G  =  —  1 +  ...,  i.e.  z  is  a  pole  of  the  Green 

(k-z)r 

function  G,  Gj^  =  0, 


G  -  g  - 
3^ 


g(x,t,l)p  (t,y,k) dt,  g  = 


exp(ikrxy) 


4Tir. 


(26) 


xy 


]j  .. 


3N4 


Multiplying  (26)  by  (k  -  z)r  and  taking  k  »  z,  we  get 


g(x,t,z) 


3R(t,y) 
9N . 


dt  =  0 , 


x  e  r 


(27) 


The  kernel  R(t,y)  is  degenerate.  Thus  a  function  a  £  0  exists 
such  that 

g(x,t;z)  0(t)dt  =  0  (28) 

•*  r 
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This  mcras' s  that  k  •••■  z  is  a  so ro  of  some  of  t ho  i'une ;  ions 

(k).  Coavorsclv,  .if  !  0  is  a  solution  of  ( 2 !) )  ,  thon 
n  ' 

u  --  [  g  ( >: ,  t. ;  z )  s‘,  ( t )  d  t  (29) 

■T 

is  a  solution  of  the  problem 

(y2  +  z2)u  =  0  in  o  ,  u  I  =  0  (30) 

1 

with  the  outgoing  asymptotic  at.  infinity.  Bonce  u  =  0  in  Q 

2 

if  7.  is  not  a  pole  of  G.  Since  z  is  complex  ,  u(x)  z  0 
in  D  (as  a  solution  of  the  homogeneous  interior  problem) . 

By  the  boundary  value  jump  relation  v  =  0.  This  contradiction 
proves  that,  z  is  a  pole  of  G.  A  variational  method  for  cal¬ 
culation  of  eigenvalues  of  nonself ad joint  compact  operators 
is  given  in  [6] . 


j 4 .  Open  problems . 

1)  The  inverse  problem  formulated  in  Section  2  is  of 
interest.  It  if;  very  interesting  to  have  partial  answers: 
what  information  about  the  geometry  of  a  scat borer  can  bo 
obtained  from  the  location  of  the  complex  polos. 

2)  There  is  a  conjecture  |3]  that  tne  coirplcfx  poles  of 
the  Green  function  of  the  problem  (i) -(4)  for  a  convex  smooth 
compact  boundary  are  simple.  It  would  bo  interesting  to  prove 
it  or  to  give  a  counte  romrnyL  • . 
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li  would  bo  interest  ino  to  perform  numerically  method 


(H'jor 


•d  in  Section  3.  b  in  some  practical  problems . 


4) 


In  [It |  none  properties  o£  the  purely  real  polos  Ke  p_ 


lm  *  0  v.'ere  established.  It.  would  be;  interesting  to  tell 
'  1 

what  information  about  the  geometry  of  an  obstacle  can  be 
obtained  fro;.:  i.  i : o  location  of  the  purely  real  poles.  In  the  ph 
literature  the  complex  plane  k  =  ip  is  usually  used.  On  this 
plane*  the  pur. -ly  real  com;  lex  poles  are  purely  imaginary, 


Re  k  . 
J 


0,  Ira  }:  .  <  0. 
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